Stray-field imaging of magnetic vortices with a single diamond spin 
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Nanostructured ferromagnetic thin films exhibiting a vortex state are of great interest both for 
fundamental studies^ and for applications such as no n- volatile magnetic storage^ ^ and microwave 
generation^. Although magnetic vortex structures have been investigated using a wide range of 
microscopy techniques^ obtaining quantitative information that is directly comparable to theory 
remains a challenging task. Further, imaging the vortex core, which can be as small as 10 nm and 
plays a crucial ro le in the vortex dynamics, is a long-standing goal that has been reached by very 
few methods onl}0^^. Here we use a recently reported scanning probe magnetometer based on a 
single Nitrogen- Vacancy (NV) defect in diamond^^ to quantitatively map the stray field above a 
thin ferromagnetic square in a vortex state. Owing to an unprecedented access to weak magnetic 
fields combined with an atomic-sized detection volume, scanning NV magnetometry reveals in 3D 
the full structure of the magnetic field distribution, including the detection of the vortex core. 
Furthermore, we show that the vectorial and quantitative nature of the measurement enable direct 
comparisons with micromagnetic simulations. This work thus opens new avenues for nanomagnetism 
and spintronics studies. 



Although a remarkable number of magnetic mi- 
croscopy techniques have been developed over the last 
decades, imaging magnetism at the nanoscale remains a 
challenging task since it requires a combination of high 
spatial resolution and sensitivit}^^^. A first approach 
consists in directly mapping the sample magnetization, 
which implies sending and collecting back test particles 
whose interaction with matter has a magnetization- 
dependent term. This approach provides the highest 
spatial resolution to date, down to the atomic-scale for 
spin-polarized scanning tunneling microscop}'^ and about 
10 nm in transmission X-ray microscopy^^. However, 
these techniques require highly complex experimental 
apparatus and a dedicated sample preparation so that 
the particles can reach and escape from the tested region 
without perturbations. To observe magnetic samples 
in their real, unprepared state, a more suited approach 
consists in mapping the magnetic stray field generated 
outside the sample, although the spatial resolution is 
then limited both by the probe size and its distance 
to the sample. Among many stray field microscopy 
techniques'^, magnetic force microscopy (MFM) has 
become ubiquitous, as it provides a spatial resolution 
below 50 nm^^ and operates under ambient conditions 
without any specific sample preparation. It was for 
instance the first method that allowed the observation 
of the core of a magnetic vortex in a thin ferromagnetic 
dislP. However, due to the intrinsic magnetic nature 
of the probe, MFM is known to be perturbative and 
not easily quantitative^', therefore limiting its field of 
applications. Recently, a magnetometer based on the 
magnetic response of a single nitrogen- vacancy (NV) de- 
fect in diamond has been proposecP^K^^, which promises 
significant advances in magnetic imaging. Indeed, it 
provides non-perturbing and quantitative measurements 
of the stray magnetic field, with an unprecedented access 



to low field s combined with an atomic-sized detection 
volumd^^. Here, we use NV-based magnetometry to 
reveal in 3D the stray field emanating from magnetic 
vortices in nanostructured ferromagnetic thin films. 
We measure several field components and demonstrate 
scanning NV magnetometry as a new method for direct 
observation of the vortex core. Reliability of the method 
is highlighted by the good agreement obtained with 
micromagnetic calculations. 

The scanning NV magnetometer combines an optical 
confocal microscope and an atomic force microscope 
(AFM), all operating under ambient conditions (see Sup- 
plementary Information). As sketched in Figure [l](a), 
a diamond nanocrystal hosting a single NV defect is 
grafted at the apex of the AFM tip and used as an 
atomic-sized magnetic senso/'^. The NV defect ground 
state is a spin triplet with a zero- field splitting D = 2.87 
GHz between a singlet state rris = and a doublet 
rris = ±1, where rris denotes the spin projection along 
the intrinsic quantization axis of the NV defect unv 
(Fig. [TJa)). Under optical illumination, the NV defect is 
efficiently polarized into the rris = spin sublevel and 
exhibits a spin-dependent photoluminescence (PL)P^. 
These combined properties enable the detection of 
electron spin resonance (ESR) by optical mean^^. A 
typical ESR spectrum of the NV sensor placed in a static 
magnetic field is shown in Figure [TJb). In the limit of 
weak magnetic fields (< 3 mT), the spin quantization 
axis remains fixed by the NV defect axis, and the two 
ESR frequencies are given hy ur = D ± 7e5Nv/27r, 
where 5nv is the magnetic field projection along the 
NV axis and 7e is the electron gyromagnetic ratio 
(Fig. [ijb)). Note that strain-induced splitting of the 
mg = ±1 spin sublevels has been omitted for clarity 
purpose (see Supplementary Information). Measurement 
of the magnetic field through Zeeman shifts of the ESR 
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Figure 1. Principle of the experiment. (a)-The apex of an 
AFM tip is functionalized with a 20-nm diamond nanocrystal 
hosting a single NV defect. A microscope objective placed on 
top of the AFM tip is used both for exciting and collecting 
the NV defect spin-dependent PL. A microwave (MW) field 
is generated by an antenna approached in the vicinity of the 
NV defect. The NV defect quantization axis unv, which gives 
the measured magnetic field component, is described by the 
angles and in a spherical coordinate system. (b)-ESR 
spectrum of the magnetic probe recorded by monitoring the 
NV defect PL intensity while sweeping the frequency of the 
MW field. The ESR splitting is proportional to the projection 
of the magnetic field along the NV axis, here Bnv = 0.9 mT. 
(c)-Typical AFM image of a thin square of Fe2oNi8o- The 
white arrows indicate the curling magnetization of the vortex 
state. 



by the NV defect is 5nv,i = ±27r(i^i — D)/^^^ and 
negative if vr = i-e. for 5nv,2 = ±27r(i/2 — ^)/7e 
(see Supplementary Information). The resulting image 
thus exhibits positive and negative signal regions corre- 
sponding to iso-magnetic field contours (5nv,17^nv,2), 
as well as zero-signal regions for any other field projec- 
tions. An example of such a dual-iso-B image recorded 
above a ferromagnetic square is presented in Figure [2](a), 
revealing a flower-shaped magnetic field distribution. 
Independent measurement of the NV defect orienta- 
tion with a calibrated magnetic field indicates that 
the NV defect is oriented mainly along the y axis in 
this experiment. Hence the magnetic map shown in 
Figure [2|a) corresponds to the field component that 
is parallel to the sample surface and along a side of 
the ferromagnetic square. By properly selecting the 
orientation of the single NV defect grafted at the apex 
of the AFM tip, any field component can be similarly 
measured. For instance, Figure [2|b) is essentially a map 
of the 2:-component, i.e. the out-of-plane magnetic field. 
These experiments establish NV-based magnetometry as 
a unique instrument providing quantitative and vectorial 
magnetic field images at the nanoscale. 

The measured stray field arises from Neel domain 
walls at the square diagonals which induce volume 
charges with opposite sign on each side of the wall. The 
general lobe structure of the magnetic distribution can 
therefore be qualitatively understood by considering 
the magnetic field created by an assembly of magnetic 
dipoles placed along the diagonals. To get more precise 



frequency is therefore quantitative and non perturbing, 
since the dipolar field from a single NV defect is as low 
as 1 /iT at 10 nm distance^^. 

In the following, we use the scanning NV magne- 
tometer to map the magnetic stray field generated by 
magnetic vortices in ferromagnetic thin films. More 
precisely, we investigate squares of Fe2oNi8o with a 
thickness of 50 nm and a 5 /im side length. This 
magnetic structure is characterized by a curling in-plane 
magnetization with a vortex core in the center where 
the magnetization points out of the plane (Fig. [ijc)). 
Such a magnetic vortex is an ideal object for evaluating 
the sensitivity and spatial resolution of any advanced 
magnetic microscopy technique since the size of the 
vortex core can reach 10 nnP^'. Furthermore this 
magnetic structure is topologically stable and therefore 
appears as an interesting candidate for memory cells in 
non-volatile data-storage deviceP^^. 

The NV probe is first approached at a distance d ^ 300 
nm above the magnetic vortex structure. Magnetic field 
mapping is performed by scanning the sample while 
measuring the difference of NV defect PL intensity for 
two fixed microwave (MW) frequencies and z^2, applied 
consecutively at each point of the scail^. This signal is 
positive if = i^i^ i.e. when the local field experienced 
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Figure 2. Stray field imaging of a square of Fe2oNi8o. 

(a),(b)-Dual-iso-B images recorded with a probe-to-sample 
distance d — 400 nm in (a) and d — 300 nm in (b). The 
red squares depict the 5 /xm magnetic structure. In (a), the 
NV is oriented mainly along the ^-axis {0 — 85°, Lp — 85°) and 
in (b) the NV defect has a strong component along the z-axis 
(6* = 37°, = 135°). Images correspond to 200 x 200 pix- 
els, with a 30-nm pixel size and a 60 ms acquisition time per 
pixel. (c),(d)-Simulated magnetic field maps corresponding 
to (a) and (b), respectively. 
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predictions, the sample magnetization distribution was 
first calculated through micromagnetic simulations using 
OOMMF softwar^, with a cell size of 5 X 5 X 5 nm'^. 
Once the equilibrium magnetization state found, the 
magnetic field generated by the structure is computed 
by summing the contribution of all magnetization 
cells. This magnetic field is then projected along the 
experimentally measured NV axis, and the NV defect 
ESR response is finally applied to get the simulated 
dual-iso-B image (see Supplementary Information for 
details). As shown in Figure |2|c)&(d), the simulations 
reproduce well the experimental data, including apparent 
asymmetries which stem from the imperfect alignment 
of the corresponding NV defects with respect to the y 
and z axis. By varying the probe-to-sample distance 
(i, it is even possible, as shown in the Supplementary 
Information, to obtain an overview of the stray field in 
the half-space above the sample, thus providing a 3D 
mapping of the magnetic field distribution. 

The magnetic images shown in Fig. [2|a)&(b) are 
recorded with a probe-to-sample distance of several 
hundreds of nanometers. At such distances, the stray 
field linked to the out-of-plane vortex core magnetization 
cannot be detected. The vortex center therefore appears 
dark (zero field) both in the images and in the simula- 
tions. With the aim of observing the vortex core, the 
NV probe is brought as close as possible to the sample 
with the AFM operating in tapping mode. Figure |3ja) 
shows the dual-iso-B image of the whole square obtained 
with the NV defect probe used in Fig. [2|b), i.e. mainly 
oriented along the z axis. By zooming in the center 
of the structure, the vortex core is revealed with an 
apparent size Lcore ^ 100 nm (Fig. |3jb)). Although 
the spatial resolution of scanning NV magnetometry is 
ultimately given by the atomic-sized detection volume, 
the effective resolving power is rather limited by the 
probe-to-sample distance, a common feature of any stray 
field microscopy technique. Indeed, numerical calcula- 
tions indicate that the size of the out-of-plane vortex 
core magnetization is around 20 nm at the sample sur- 
face. The probe-to-sample distance can be estimated by 
comparing the experimental results to magnetic images 
simulated at different distances from the sample surface. 
As shown in Figures [3]^c),(d)&(e), a good agreement is 
obtained for a distance d ~ 100 nm. We attribute this 
relatively large value to an imperfect positioning of the 
diamond nanocrystal at the apex of the AFM tip. A 
more precise control of the NV magnetic sensor position 
could be achieved by using diamond nanopillar-probes^^. 
Although the simulation fairly agrees with experimental 
data, we note that the sharp structures lying along 
the square diagonals (Fig. [sj^cV) are not observed in 
the experimental images (Fig. |3[a)). In those regions, 
the stray field generated by Neel domain walls is 
strong enough to induce a mixing of the NV defect 
electron spin sublevels which results in an overall reduc- 
tion of ESR contrast^^ (see Supplementary Information). 
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Figure 3. Observation of the vortex core. (a)-Dual-iso-B 
image of the whole square obtained with the AFM operat- 
ing in tapping mode. The same NV probe as in Fig. is 
used {0 = 37°, (p — 135°). (b)-Zoom into the center of the 
structure enabling the detection of the stray field emanating 
from the vortex core magnetization. The image corresponds 
to 110x110 pixels, with a 4.7 nm pixel size and a 60 ms acqui- 
sition time per pixel. (c),(d)-Simulated magnetic field maps 
corresponding to (a) and (b), respectively, with a probe-to- 
sample distance d — 100 nm. (e)-Line-cut taken along the 
red dashed line in (b) (markers) and (d) (solid line). (f),(g)- 
Dual-iso-B image of the square center while applying a bias 
field with a projection along the NV axis Bext = — 0.8 mT in 
(f) and Bext = +0.8 mT in (g). 



The vortex structure is commonly characterized 
by two independent binary properties. The chirality 
c determines whether the in-plane magnetization is 
curling clockwise (c = +1) or counterclockwise (c = — 1), 
while the polarity p of the vortex indicates the upward 
{p = +1) or downward {p = —1) magnetization of the 
core. The apparent pairing of the zero-field contours 
around the vortex core observed in Figure IsFb) directly 
enables us to determine that c x p = +1 (see Supple- 
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Figure 4. Higher-order magnetization distribution. 

(a)-Magnetic field distribution recorded above a square of 
Fe2oNi8o for a probe-to-sample distance d = 300 nm. The NV 
defect orientation is given by ^ = 65° and cp = 30°. (b)-AFM 
image of the square of Fe2oNi8o- The white arrows depict 
the higher-energy magnetization structure, with four mag- 
netic vortices (red dots) distributed around one anti-vortex 
at the center (blue dot). (c)-Simulated image considering a 
magnetization following the structure shown in (b). 

mentary Information). We note that only an absolute 
value of the field is measured, preventing a direct 
measurement of the core polarity without additional 
measurements. This limitation could be overcome by 
addressing selectively one of the two ESR transitions 
through circularly polarized MW excitatioii^. Here, 
we gain further insights into the vortex structure by 
applying an external bias magnetic field along the 
z-axis with a projection ^ext along the NV axis. The 
amplitude of this external magnetic field is chosen weak 
enough to avoid modifying the sample magnetization, 
so that 5ext just adds up to the regular stray field of 
the vortex structure. Figures [3]^f)&(g) show magnetic 
images obtained with ^ext = —0.8 mT and ^ext = +0.8 
mT, respectively. The stray field of the vortex core 
turns dark - the net magnetic field vanishes - with the 
negative bias field, which reveals that the vortex core 
magnetization points upward p = +1, and indicates a 
stray field emanating from the core ^core ^ +0.8 mT. 



Conversely, the core field is increased with a positive 
bias field since it adds up to ^core- From this polarity 
measurement, the chirality c = +1 can finally be 
deduced leading to a full characterization of the vortex 
structure. 

As a final experiment, we study a higher-order 
magnetization distribution on a Fe2oNi8o square dot. 
As shown in Figure [4]^a), highly non trivial magnetic 
field distributions are observed above some square 
ferromagnetic structures of the sample. Obviously 
these particular squares are not in the regular vortex 
distribution studied previously. One way to determine 
the magnetization distribution is to make an assump- 
tion about it, compute the expected dual-iso-B-image, 
and compare to experimental data. Here the four 
near-symmetry points distributed around the square 
center suggest a magnetization state comprising four 
vortices distributed around one anti- vortex at the center 
(Fig. gb)), which is a known equilibrium state in such 
ferromagnetic structures'^. Comparison between the 
simulated (Fig. ^c)) and measured magnetic images 
unambiguously validates our hypothesis that this par- 
ticular square exhibits a four-vortex-one-antivortex 
magnetization state'^. 

This experiment illustrates that scanning NV magne- 
tometry enables robust comparisons with micromagnetic 
simulations and therefore appears as a powerful tool 
for fundamental studies in nanomagnetism and, more 
generally, nanoscience. Indeed, the direct measurement 
of weak magnetic fields, without perturbation and with a 
nanoscale resolution, may allow clearing some important 
issues, such as the nature of magnetic domain walls in 
ultrathin films with perpendicular anisotropy, which 
controls their current-induced motion'^. The detec- 
tion of exotic magnetic structures known as skyrmion 
lattices^^ is another example of a challenge that may 
be met by scanning NV magnetometry. On a broader 
perspective, the observation of the stray fields created by 
persistent currents in micro- and nanostructures where 
the phase coherence of the carriers is preserved, a quan- 
tum realization of the naive Amperian currents proposed 
two centuries ago to explain the magnetism of matter, 
appears also within reach of scanning NV magnetometry. 
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SUPPLEMENTARY INFORMATION 
A. Experimental setup and sample fabrication 

The experimental setup combines a tuning-fork-based atomic force microscope (AFM) and a confocal optical mi- 
croscope (attoAFM/CFM, Attocube Systems). A detailed description of the setup as well as the method to graft a 
diamond nanocrystal at the apex of the AFM tip can be found in Ref.[lJ. We use commercially available diamond 
nanocrystals (SYP 0.05, Van Moppes SA, Geneva), in which single NV defects were created through high energy 
electron irradiation (13.6 MeV) followed by thermal annealing at 800° C. The irradiated nanocrystals were finally 
oxidized in air at 550° C during two hours. This procedure enables to reduce the size of the nanodiamond^ and 
leads to an efficient charge state conversion of the created NV defects into the negatively-charged statJ^. For the 
experiments reported in the main article, the size of the nanodiamonds were in the ^ 20 nm range, as verified using 
AFM images before grafting the nanodiamond at the apex of the tip. The unicity of the NV defect was checked 
through measurements of antibunching in the second order correlation function of the NV defect photoluminescence 
(PL), using a standard Hanbury Brown and Twiss interferometei^. 

For electron spin resonance (ESR) spectroscopy, a microwave (MW) excitation is applied through a 20 fim copper 
wire directly spanned on the magnetic sample. The 5x5 /im^ square dots were patterned on a silicon substrate using 
electron-beam lithography followed by evaporation of 50 nm of Fe2oNi8o and lift-off. All experiments were performed 
at room temperature. 



B. ESR frequencies vs. magnetic field 

In this section, we give the analytical expression of the NV defect ESR frequencies as a function of the applied 
magnetic field, while taking into account strain effects of the diamond lattice. 
The ground- state spin Hamiltonian of the NV defect electron spin (5 = 1) reads 

n = hDSl + hE{Sj^ - S^) + n-feB • S (1) 

where h is the Planck constant, D is the zero-field splitting, E the strain splitting, Z the NV defect axis, 7e the 
electron gyromagnetic ratio, and B the local magnetic field [Fig. [sja)]. By numerically computing the eigenenergies 
of the two ESR frequencies iyR^± can be calculated for any magnetic field B. For instance, the ESR frequencies 
^R,± are plotted in Figure [sjb) as a function of the magnetic field amplitude B = ||B|| for various angles a between 
B and the NV defect axis. 

With the aim of obtaining a simple analytical expression for vr^±^ the Hamiltonian 1-L is written 

n = Hll+H±+ ^strain (2) 

with 

7^11 =hDSl + h-ieB^wSz 

ni. =h^e{BxSx + BYSY) (3) 
'^strain = hE{Sj^ — Sy) , 

where we introduced the notation B^y = Bz as used in the main article. Typical values for the zero-field parameters 
in diamond nanocrystals are D ^ 2.87 GHz and E ^ 5 MHz. One can therefore distinguish different regimes depending 
on the value of B\\ = \B^y\ and B_i = B\ + By , which are the parallel and transverse components of the magnetic 
field with respect to the NV axis, respectively. 

Focusing on weak magnetic field amplitudes, we consider that <C Hy. The ESR frequencies are then simply 
given by 

z^i?,±(5Nv) =D± V(7e5Nv/2^)2 + E^ . (4) 

This formula can be used whenever the transverse component of the magnetic field B^ is much smaller than 2i\D j^^ ^ 
100 mT. This is illustrated in Figure |5|b), where the prediction of Equation Q is plotted in dotted lines together 
with the full calculation (solid lines). Obviously the approximate formula is relevant for field amplitudes in the range 
0-2 mT, which is the range we address in this work. In addition, if jeB\\/27r E, then ^strain can be neglected and 
the ESR frequencies become 



^i?,±(^Nv) = D± 7e5Nv/27r , 



(5) 
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Figure 5. (a)-Schematic of the NV defect. A magnetic field B is applied with an angle a with respect to the NV defect axis. 
(b)-ESR frequencies iyR,± as a function of the magnetic field amplitude B = ||B|| for different angles a. The solid lines are 
obtained by diagonalizing the full Hamiltonian described by Equation ([T]), while using D = 2.87 GHz and E — b MHz. The 
dashed lines are obtained by using Equation Q. (c)-Optically-detected ESR spectra of a single NV defect placed at the apex of 
the AFM tip for Bnv = mT and Bnv = 1.2 mT (down-shifted for clarity). The solid lines are data fitting with Equation ([6| 
which allows us to determine the ESR linewidth Az/ =13 MHz and the ESR contrast C — 0.11, for this particular NV defect. 
(d)-ESR frequencies as a function of Bnv, with a — 44°. The solid lines are the predictions of Equation Q with the zero- field 
parameters D = 2869.5 MHz and E = 10.5 MHz. 



as indicated in the main paper. 

Experimentally, the ESR frequencies are inferred from the optically-detected ESR spectrum, i.e. the NV defect PL 
rate IZiv^ ^nv) as a function of the MW frequency v. Once normalized to unity, it can be expressed as 



7^(^/,5Nv) = l-C 



Ai//2Vlri2 



Ai//2vlri2 



(6) 



where G{x) = e~^ denotes a Gaussian function, /S.v is the ESR linewidth (FWHM) and C its contrast. Two measured 
ESR spectra are shown in Figure |5|c) for a single NV defect attached to the AFM tip, together with data fitting 
using Equation ([6|. Varying 5nv enables us to verify Equation Q experimentally, as illustrated in Figure [5|d). Note 
that the vanishing slope at very low magnetic fields - i.e. when ^^3/271 <^ E - is responsible for the fact that in our 
dual-iso-B images, the iso-B lines corresponding to zero field are always wider that those corresponding to a non-zero 
field. 



C. Determining the orientation of the NV defect 



The orientation unv of the NV defect gives the projection axis of the measured magnetic field in all magnetometry 
experiments, i.e. 5nv = B • unv- To determine the spherical angles and (f that characterize unv [Fig. [oj^a)], we 
apply a calibrated magnetic field Bcoii = ^coiiUcoii using a three-axis coil system [Fig. ^h)] while monitoring one 
of the two ESR frequencies of the NV defect {e.g. The procedure we employed is as follows. We first set 

^coii = 2 mT and ^coii = 90°, and sweep the angle (/?coii from to 90°. This yields the angle cp by fitting Equation 
Q to the data, using the fact that 5nv cos((/?coii — ^) [Fig. [sjc)]. We next set (fcoH = ^ and sweep ^coii from 
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to 90°. The angle 6 is deduced by fitting Equation Q to the data, given that 5nv ^ cos(6>coii — 0) [Fig. |6jd)]. This 
procedure enables us to determine the angles and in the range [0, tt], which completely defines the direction of unv- 



a ^ 



c 

r 2.92 



X 



2.5 



a: 

C/D 
LU 




b ^ 







|(/p = 70±0.8°| V 



90 

(/^coii (deg.) 



180 




N 
X 


2.92 






>. 




o 




CD 






2.90 


D- 
0) 








C/5 




LU 


2.88 









V |6' - 120±0.2°| 



90 

6'coii (deg.) 



180 



Figure 6. (a)-The axis of the NV defect unv is described by the angles and Lp in the spherical coordinate system of the 
laboratory reference frame. (b)-An external magnetic field is applied using a three-axis coil system. Its direction, of unit vector 
Ucoii, is described by the angles ^coii and (/?coii in the same coordinate system. (c)-Measured ESR frequency fR^^ as a function 
of the angle ^coH of the applied magnetic field, while fixing ^coii = 90°. The solid line is a fit to the data using Equation Q, 
which yields the angle Lp. (d)-Measured ESR frequency z/h,+ as a function of the angle ^coii while keeping LpcoH = ^- The solid 
line is a fit to the data using Equation Q, which yields the angle 0. 



D. Magnetization and stray field simulations 



The magnetization distribution inside the square film of Fe2oNi8o was computed using OOMMF softwar^. The 
parameters used for Fe2oNi8o were: spontaneous magnetization Mg = 810^A/m and an exchange constant A = 
10~^^ J/m. A magnetization cell size 5x5x5 nm^ was used, which is below the micromagnetic exchange length in 
order to describe the vortex core structure precisely. From the calculated equilibrium magnetization [Fig. [7|^a)], the 
stray magnetic field is computed by summing the contribution of all magnetization cells, and then projected along the 
NV defect axis in order to reach a map of 5nv- As an example. Figure [Tj^b) shows the calculated z-axis component 
of the stray field for a probe-to-sample distance d = 200 nm. 

In the following, we take into account the strain-induced splitting, which was omitted in the main text for clarity 
purpose. We therefore use Equation Q for the simulation of the magnetic field maps recorded with the scanning 
NV magnetometer. The dual-iso-B images shown in the main article are obtained by measuring the difference of NV 
defect PL intensity for two fixed MW frequencies ui and U2- The resulting signal reads 



5(i^i,i^2,^Nv) =7^(^^2,^Nv) -7^(^^l,5Nv) 



(7) 



where 7^(i/, 5nv) is given by Equation ([6|. The signal S is positive (resp. negative) when vi (resp. 1^2) 
matches one of the two ESR frequencies iyR^± - i.e. when B^y = 5nv,i = ±27rY^ {ui — DY — E'^ /^^ (resp. 

^NV,2 ±27l^f{h^ 



DY — E'^ Me) ~ and is null otherwise. 
As an illustration, we consider the same NV defect as in Figures [5|c)&(d), with an ESR spectrum at zero magnetic 
field reproduced in Figure [tJ^c). We choose vi = 2905 MHz and 1^2 = 2880 MHz, which corresponds to 5nv,i = 1-2 
mT and 5nv,2 = 0. Having fully characterized the ESR response, i.e. its linewidth and contrast, the resulting 
dual-iso-B signal tS(i^i, i^2, ^nv) can be computed as a function of ^nv, as shown in Figure 7|^d). The simulated 
magnetic field maps are then obtained by applying the i^2, ^Nv) response to the calculated stray field once 

projected along the NV defect axis. An example of a calculated dual-iso-B image is shown in Fig. [7|^e), where we 
assumed d = 200 nm and the NV axis to be along the z axis, i.e. Z = z. 



9 




Figure 7. (a)-Map of the in-plane magnetization of a 5 /xm Fe2oNi8o square dot, as computed by OOMMF software, (b)- 
Calculated distribution of the Bz component of the stray magnetic field at a distance d = 200 nm above the sample surface. 
(c)-ESR spectrum measured at Bnv = 0. The microwave frequencies applied for the dual-iso-B images are chosen to be 
= 2905 MHz and 1^2 = 2880 MHz, as indicated with blue arrows. (d)-PL difference signal S {ui , 1^2 , B-^y) as a function of 
Bnv, calculated from the ESR spectrum in (c) and Equation ([7|. (e)-Dual-iso-B image calculated by applying S {vi , 1^2 , B-^y) 
to the magnetic field map shown in (b). Here we assumed the NV axis to be parallel to the z axis. 

E. Distance dependence 

Magnetic field images recorded above a 5 /im Fe2oNi8o square dot at several probe-to-sample distances, varying 
between d = 100 and d = 600 nm, are shown in Figure [Sj For these experiments, the probe-to-sample distance is 
controlled by using the AFM in a dual-pass mode with a constant lift height during the second pass. The same NV 
defect probe as in Figure 3 of the main article is used, i.e. mainly oriented along the z-axis. A good agreement 
with the simulations is found, which illustrates the three-dimensional stray field mapping capabilities of scanning 
NV magnetometry. The slight asymmetries observed in both the experimental and simulated images result from the 
imperfect alignment of the NV defect axis with respect to the z-axis. 



100 nm 250 nm 400 nm 600 nm d 

\ \ \ \ ^ 




5 /im 



Figure 8. Experimental dual-iso-B images (top panels) obtained for a NV defect oriented mainly along the z axis {0 = 37°, 
(p = 135°) at different probe-sample distances together with the simulations (bottom panels). The red square depicts the 
position of the 5 /xm square dot. 
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F. Magnetic-field-induced fiuorescence quenching 



Although the simulations fairly agree with experimental data, the sharp structures lying along the square diagonals 
are not observed in the experimental image recorded at the smallest probe-to-sample distance d — 100 nm [Fig. [sj. 
As explained in the main text, a strong stray field is generated by Neel domain walls in these regions. Indeed, the 
calculated magnetic field maps at d = 100 nm indicate that the stray field components along the y- and z-axis reach 
nearly ~ 10 mT at the square diagonals, as illustrated in Figures |9ja)&(b). In this condition, the approximation 
V.^ ^11 is not fulfilled [see Eq. ([2|] and the quantization axis is determined by the local magnetic field rather than 
by the NV defect axis. As a result, the spin projection along the NV axis is no longer a good quantum number 
and the eigenstates of the spin Hamiltonian are given by superpositions of the = and = ±1 spin sublevels''. 
In this regime, the contrast of optically detected ESR vanishes^. Magnetic field imaging through optically-detected 
ESR is therefore inefficient in the regime of 'high' off-axis magnetic field. This explains why the zero-field lines are 
not visible at the square's diagonals for d = 100 nm, where a strong off-axis magnetic field is experienced by the NV 
defect. 

It was also shown that the decreased ESR contrast is accompanied by a reduction of the NV defect PL intensity 
when the magnitude of increases^ ^. Such a magnetic-field-dependent PL can be exploited as an alternative 
method to map high off-axis magnetic field region^^. This is illustrated in Figure [9|c) where dark lines are observed 
along the square's diagonals in the PL image for a probe-to-sample distance d = 100 nm. In this experiment, we 
record the NV defect PL intensity without applying any MW field. The dark lines indicate a PL quenching induced 
by the off- axis component of the magnetic field generated by Neel domain walls. Near the center of the magnetic 
structure, the magnetic field amplitude decreases [see Fig. [9]^c)], which enabled us to image the vortex core using the 
ESR response of the NV defect. 

We note that the resulting PL quenching image exhibits obvious similarities with the one recorded with a magnetic 
force microscope (MFM), which is sensitive to the field gradient [Fig. loFd)]. 



By (mT) (mT) PL (a.u.) 




MFM image 



Figure 9. (a),(b)-Calculated By and Bz components of the stray field distribution at a distance d — 100 nm. (c)-PL quenching 
image obtained without applying any MW field and with the AFM operating in tapping mode. Zooming into the center shows 
that the PL quenching vanishes in the vicinity of the center of the structure, where the magnetic field amplitude is expected 
to decrease. (d)-MFM image of a similar square of Fe2oNi8o- 



G. Vortex polarity and chirality 



The product of the vortex chirality c and the vortex core polarity p can be determined by simply examining the 
zero-field contours around the vortex core in the experimental dual-iso-B images. Indeed, zero-field lines separate the 
regions of positive and negative magnetic field projection. It is then possible to determine which "lobes" have the 
same sign as the vortex core field. Figure [lOja) shows the magnetic field map calculated at d = 100 nm from the 
sample surface and projected along the z-axis, for all possible combinations of c and p. As expected, the pairing of 
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the zero-field contours around the vortex core (black dotted lines) are fixed by the product c x p. Comparing the 
experimental results to simulations directly enables to determine that c xp = +1 for the vortex state studied in Figure 
3 of the main text [reproduced in Fig. [lOjb)]. 

By applying an external magnetic field B^^t that cancels out the net field above the vortex core, the polarity of 
the vortex core can be determined as shown in the main paper. In Figures [Tojc), we rather show dual-iso-B images 
obtained for a bias field ^ext antiparallel to the vortex core's field ^core such that |5core| — l^extl is slightly positive 
or slightly negative. In the latter case, the net field above the vortex core is now of opposite sign compared to the 
zero- field case, which reverses the pairing of the zero- field lines in the image. 



c = +1 ,p = +1 c=-l,p = -hl c X p = +1 




Figure 10. (a)-Calculated Bz component of the stray field produced by the center of the magnetic structure at d = 100 nm, 
for the four possible combinations of chirality c and polarity p of the vortex state. The black dotted lines indicate the zero- 
field contours, separating positive (upward) and negative (downward) field regions. (b)-Magnetic field map of the stray field 
emanating from the vortex core magnetization, as shown in Figure 3 of the main text. (c)-Dual-iso-B image of the square 
center while applying a bias field with a projection along the NV axis Be^t = — 0.6 mT and Be^t = — 1.3 mT. 
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